Pion Number Fluctuations and Correlations 
in the Statistical System with Fixed Isospin 

V. V. Begun, 1 ' 2 M. I. Gorenstein, 1,2 and 0. A. Mogilevsky 1 

1 Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 
2 Frankfurt Institute for Advanced Studies, Frankfurt, Germany 

Abstract 

The statistical system of pions with zero total isospin is studied. The suppression effects for the 
average yields due to isospin conservation are the same for ir°, ir + and ir~ . However, a behavior of the 
corresponding particle number fluctuations are different. For neutral pions there is the enhancement 
of the fluctuations, whereas for charged pions the isospin conservation suppresses fluctuations. The 
correlations between the numbers of charged and neutral pions are observed for finite systems. This 
causes a maximum of the total pion number fluctuations for small systems. The thermodynamic limit 
values for the scaled variances of neutral and charged pions are calculated. The enhancements of the 
fluctuations due to Bose statistics are found and discussed. 
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I. INTRODUCTION 



The aim of the present paper is to study some aspects of non-Abelian symmetries in the 
statistical models. We consider SU(2)-isospin symmetry group for the pion system. The role of 
the isospin conservation in a many body system was first considered in the pioneering paper of 
Bethe [1|. Many efforts were then aimed at studies of the pion system with fixed isospin 
An effective theoretical formalism for non-Abelian symmetries in the statistical mechanics was 



developed in Ref. J4] on the basis of the group projection technique. It allowed to consider the 
impact of the isospin conservation on the particle abundances and the form of their momentum 
spectra in the statistical models of hadron production js , 6 ] . The group projection technique was 
also used to calculate the colorless partition function of the quark-gluon gas with SU(N c )-color 
symmetry [3, Is]. 

Our primary interest in the present paper is to study an influence of 
conservation on the particle number fluctuations. It was recently found 
vation of Abelian (additive) charges causes the suppression of the particle number fluctuations. 
In the present study we restrict our consideration to the simplest statistical system with non- 
Abelian symmetry - an ideal pion gas with zero isospin 1 = 0. Most discussions are done 
within Boltzmann statistics. This makes possible to obtain transparent analytical results and 
compare them with those in the canonical ensemble for zero electric charge Q = 0. 

The paper is organized as the following. In Section [III we consider the partition function and 
total number of pions (the mean value and scaled variance). In Section HTT1 we calculate the 
mean multiplicities, fluctuations and correlations for 7r + , 7r~ and 7r° mesons. In Section HVl we 
consider the specific effects due to Bose statistics. The Section [V] summarizes the paper. 

II. PARTITION FUNCTION AND TOTAL NUMBER OF PIONS 

The partition function of the ideal Boltzmann gas of pions 7r + , 7r~, 7r° in the grand canonical 
ensemble (GCE) reads, 

7 ^ (Ao zf° (A+ zf+ (A. zf- f . . 4n , m 

Z 9 ce = ^ — ^ JT\ am = exp[(A + A + + A_)z] , (1) 

N ,N+,N-=0 u ' +' 



non-Abelian charge 
9] that exact conser- 
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where z is the one-particle partition function, 

V f°° a ( vV + m 2 \ V 



2tt 5 



Here and T are the system volume and temperature, m is the pion mass (we neglect a small 
difference between the masses of charged and neutral pions), and K2 is the modified Hankel 
function. The auxiliary parameters Xj with j = 0, +, — are introduced to calculate the mean 
pion multiplicities, fluctuations and correlations. We take Xj = 1 in the final formulae. 

In the case of exact charge conservation, i.e. in the canonical ensemble (CE) with zero charge 
Q = 0, the partition function is (see, e.g. 

Zq =° = 2. S{N+ - N ~ ] -Ivo! ivy KT" 

N ,N+,N-=0 u + 

1 /" 27r 

= exp (Ao z) — / d(j) exp [ z (A + exp[i0] + A_ exp[— i<f>]) ] . (3) 
2tt Jo 

Note that an exact charge conservation in the CE ([3]) does not affect the neutral pions. Their 
number distribution remains the Poissonian one, the same as in the GCE (JT|). 

The partition function with total isospin 1 = can be obtained using group projection 
technique. Pions are transformed under vector (adjoint) representation of the SU(2) group. 
This group has three parameters which can be chosen as Euler angles a = a, (3, 7. In this case 
the diagonal matrix elements have the following form 111 ] : 

£>i li±1 (a,/3, 7 ) = ( 1 + C 2 ° S(/3) ) , £>io(«,/3,T) = «*(J3) . (4) 

The partition function is then presented as jf]. [12I]: 

[X zD^(a)] N ° [X + z DKa)}^ [\_ z Dl^ffl* 



Z I=0 = / dfi ^2 



(5) 



iV ! N + \ NJ. 

N ,N+,N-=0 u + 

= J dfi exp [Ao z £>o )0 (<3) + A + z D\ l {a) + X^ z D^ 1 _ 1 (a 

Substituting explicit expressions for the Haar group measure dfi and matrix elements D\ 3 1 <^ 
in Eq. (jSJ), one obtains: 

j rife rife pit 

Zi=q = -—7 / da d'y d(3 sin j3 (6) 



■57T 2 



x exp 



1 _|_ cos ft s 

A -2cos/3 + z ( ) ( A + exp[i(a + 7)] + A_ exp[— i(a + 7)] 
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The change of variables 



Zi=o — — 



Arc 



a + 7, (p — (a — t)/2, cos(/3) = x and integration over gives 

1 + x 

(X + exp[i(f)} + A_exp[— 



p2tt 






/ d4 


J dx exp 


A zx + z 


/o 







(7) 



Comparing Zj = q ([7]) with the partition function Zq = q Q one observes an additional re- 
integration in Eq. (|7|). It reflects a presence of the particle number correlations between neutral 
and charged pions which were absent in the GCE and CE. 

The partition functions (JTJ) , (JH]) , and (|7|) can be simplified by taking \j = A. One finds: 

Z gce = exp(3Az) , Z Q=0 = exp(Xz) I (2Xz) , Z I=0 = exp(Xz) [I (2Xz) - h(2Xz)) , (8) 

where I n are the modified Bessel functions. The final expressions for the partition functions 
correspond to A = 1. Taking the derivatives of Z over T and V one finds the thermodynamical 
functions of the pion system. The derivatives over A give the moments of total pion number 
distribution, 



1 d 



dZ\ 



(9) 



A=l 



— v <9a; 

With Eq. one calculates the average pion multiplicity (N) and the second moment (iV 2 ) 
using the partition functions (jSJ) of different statistical ensembles: 



(N) gce 
(N} Q =o 

(N) I=0 



3z , 

I x (2z) 
■■ z + 2z ^4 • 
h(2z) ' 

h{2z) - I 2 (2z) 



(N 2 ) gce 
(iV 2 ) Q =c 



3z + 9z 2 , 

= z + 5z 2 + 4z 2ll ^ Z 



h{2z)-h{2z) ' 
The ratio Rn and the scaled variance 

(iV> 



(N 2 ) 



1=0 



(N) 



1=0 



+ z 



h{2z) ' 
2 I {2z) - h{2z) 
I Q {2z) - h{2z) 



R 



N 



X - (NY 



(10) 

(11) 

(12) 
(13) 



3z 1 (N) 
for the total number of pions calculated in different statistical ensembles are shown in Fig. [TJ 
In the GCE one obtains R N (GCE) = u N {GCE) = 1. 

Figure [T] shows that conservation laws suppress the mean particle number. For the same 
volume V and temperature T the average number of pions in the CE with Q = is smaller than 
the GCE value 3z, thus R^ < 1. Even stronger suppression effect is observed in the ensemble 
with 1 = 0, i.e. R T N < R^ < Rn{GCE) = 1. The opposite effect is seen for the scaled variance 
of the pion number fluctuations: u ! N > > un(GCE) = 1. To find the asymptotic behavior 
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FIG. 1: The ratio Rm (left) and the scaled variance ojn (right) in the GCE (the dotted lines), CE with 
Q=0 (the dashed lines), and the statistical ensemble with 1 = (the solid lines) 

of Rn and u^- at large z in the statistical ensembles with Q = and 1 = one can use the 
expansion of the modified Bessel functions at z 1 13] : 

exp(2z) 



Wz) 



TTZ 



An 2 -I / 1 

1 — + O - 



(14) 



16 z \z" 

in Eqs. ( fl0iri2l) . This gives, Rn — > 1 and % — > 1 in the thermodynamic limit, z — >■ oo, for both 
Q = and 7 = ensembles. It means that the suppression of the pion average multiplicity and 
the enhancement of the multiplicity fluctuations are the finite volume effects. 
The partition functions (0Q), (EJ), and (J7J) can be presented as, 

oo N n oo 

jV=0 7V =0 ^=0 

The degeneracy factor g(N) for different statistical ensembles can be found from Eq. (151) . To 
calculate F(iV, jV ) one needs Eqs. (JU El El). The values of g(N) and F(iV, N ) in different 
statistical ensembles are presented in Table I for N < 4, see also [2|. Note that g(N) = 3 N 
in the GCE, and g(0) = 1 in all ensembles. Table His helpful when one considers very small 
systems with only a few pions. For z < 1 one finds from Eq. (JT5J) and Table [U 



^-2 _,3 -,4 

1 + 1-. + 3.- + 7.- + 19.- 



_,3 o , ^ L 

r /= o-l + l- 2J+1- 3J+3- - . (16) 
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TABLE I: The degeneracy factors g(N) and F(N,N ) (USD in the GCE, in the CE with Q = 0, and 
in the statistical ensemble with 1 = 0. 



The mean multiplicity and higher moments can be calculated as: 



N=0 



N=0 



This gives: 



(iV) Q=0 = z + 2z 2 - z 4 



(N 2 ) 



(N) I=0 



= z 2 + 



z 

y 



Q=0 



(iV 2 )/=o = 



z + 5z 2 + 4z 3 , 



2z 2 + -z 3 + z 4 
2 



(17) 

(18) 
(19) 



The results (I18H19P can be also obtained from Eqs. (11011121) by expanding the modified Bessel 
functions at z <ti 1 [131 ]. 

r n+2 



Wz) 



+ 



n\ (n + 1)! 

For R N and u;^ f|T3|) one finds from Eqs. fjl8til9|) : 

R Q a 1 i 2 y 
k n — 3 3 ' 



O {z n+4 ) 



ryl ~ 
K N = 



3 ' 



N — 

I £J 

N — 



1 + 2z 



Az 1 



z z" 
2+ 2 - 4" 



(20) 

(21) 
(22) 



The behavior of Rn (12ip and u;^ ( |22|) at small z can be seen in Fig. [H The scaled variance ujn 
has a maximum at z « 0.5 in both Q = and 7 = statistical ensembles. 



III. FLUCTUATIONS AND CORRELATIONS OF vr , vr+, tt" 



To calculate the mean multiplicities, correlations, and fluctuations for neutral and charged 
pions one has to return back to presentations of the partition functions by Eqs. dH O [T]) . One 
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finds, 



1 dZ 



Z d\j 



(N t N 3 ) 



A=l 



1 d f dZ 



(23) 



A=l 



Using Eq. (I2"3"j) one obtains for the mean multiplicities of neutral and charged particles: 

h(2z) 



(N ) g ce 

(iV )/=o 



gee 



1=0 



' , (N ) Q=0 
z h(2z) - I 2 (2z) 



z , {N±}q =0 = z 



3 I {2z) - h{2z) 
where (N±) = (N+) = (JV_). The ratios 

(^o) D (N ± ) 



(24) 
(25) 



R 



± 



are shown in Fig. [2] (left) for Q = and 7 = statistical ensembles. 
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FIG. 2: The ratios Rq and R± (left) and the scaled variances ojq and oj± (right) in the CE with Q = 
(the dashed lines) and in the ensemble with 1 = (the solid lines). 

Note that Ro = R± = 1 in the GCE. In the CE, R® = 1, i.e. the mean number of neutral 
pions is not affected by the c harg e conservation law. The mean number of charged pions is 
suppressed in the CE, R± < 1 lOj. The behavior in the statistical ensemble with 1 = differs 
from that in the CE with Q = 0. The pion mean numbers are the same for all charge pion 
states 7T°, 7r + , and tt~ . The suppression of these pion multiplicities at I = is stronger than 
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that for 7T or 7r + in the CE with Q = 0: 

Rl = Ri < R2 < R$ = 1 . (27) 

The asymptotic behavior of (EU [25]) at large z can be found using Eq. ([TJJ. One obtains 
Rj — > 1 at z — > oo for all j = 0, +, — in both Q = and 7 = statistical ensembles. Thus, the 
suppression of the mean numbers of charged and neutral pions is the finite volume effect. At 
small z one can use Eq. (I2"0"j) . this gives at z <C 1: 

R2 - * , *J = J4*f. (28) 

The second derivatives in Eq. ( )23|) can be calculated using the partition functions (DQ El [7]) 
for different statistical ensembles. They give the second moments and correlations in the GCE: 

(N 2 ) gce = (Nl) gce = z + z 2 , (N N ± ) gce = z 2 , (29) 

in the CE: 

(N 2 ) Q=0 = z + z 2 , (N 2 ± ) Q=0 = z 2 , (N N ± ) Q=0 = z (N ± ) Q=0 , (30) 

and in the pion system with 1 = 0: 

/N 2, /at \ 1 z(7z + 4) Zq(2z) - (7z 2 + 2z + 4) A(2z) 

(iV )/=o - (iV )/=o + ^ J (2.) - ^(2.) ' (31) 

(N ± ) I=0 - (N ± ) I=0 + - io{2z) _ h{2z) , (32) 

{N ° N±)l =° ~ 15 7 (2*) - ^(2*) ' (33) 

where (N£) = (N 2 ) = (N 2 ) = (iV+AL). 

Using the above expressions one can easily construct the scaled variances Uj and correlation 
coefficients p^, 

that define the main characteristics of the pion multiplicity distributions. 

From Eqs. (I24H251 13~TH3"2"1) and asymptotic expansion of the modified Bessel functions ( 1T4"|) 
one finds the behavior of the scaled variances Uj ( l34"j) in the thermodynamic limit z — > oo. As 
seen from Fig. [2j the isospin conservation with 1 = gives the same pion number fluctuations 
as the CE with Q = 0, uJq — > 1 and u± — > 1/2, at z — > oo. 
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However, the results for finite systems are rather different. In the CE with Q = the 
fluctuations of neutral particles are the same as in the GCE, Uq = 1. The scaled variance 
u± in the CE with Q = was calculated in Ref. 9]. At z < 1 using Eq. fl20|) . one finds: 
oo2 = 1 - z 2 /2 . 

In the statistical ensemble with 1 = the scaled variances Uq and u± can be calculated for 
small system with z < 1 using Eqs. ( [2411251 13TH321) and asymptotic expansion of the modified 
Bessel functions ( 1201 . However, it is more instructive to obtain approximate expressions for 
the distributions of tt° and 7r =l= numbers using Eq. (|T5"|) and F(N, N ) values from Table [B Let 
us consider the pion states with total number of pions N < 4, i.e. we neglect small terms of 
the order of z 5 . The probability distribution P (N ) then reads: 

i 4 N 

Po(No) = yJ2 F ( N ' N °)m- 



N=0 



(35) 



For the statistical ensemble with 1 = one finds from Eq. (135]) and Table [J 
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Similar expressions for P±(N±) in the statistical ensemble with 1 = are equal to: 
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(36) 
(37) 
(38) 

(39) 
(40) 



The distributions Pq and P± are rather different. The one and two particle states have evidently 
different probabilities for neutral and (negative) positive pions. At z 1 the main configura- 
tions consist of 1 (negative) positive pion and 2 neutral pions. These states correspond to the 
total number of pions N = 2, and their probabilities are proportional to z 2 . The states with 
only 1 neutral pion have smaller probability, z 3 , as they can only appear for > 3. Inspite 
of differnces in the Pq and P± distributions they give, however, the same average numbers of 
neutral and (negative) positive pions: 



(Ar > J=0 = l-Po(l) + 2-P (2) + 4-P (4) 
(N ± ) I=0 = 1-P ± (1) + 2-P±(2) 



z 2 z 3 
= — + — 
3 6 



z 2 z 3 
= — + — 
3 6 



(41) 
(42) 



9 



Equations (141f42p clearly demonstrate that the same average numbers of 7r° and tt^ come from 
different pion states. The results for higher moments of Pq and P± distributions are different: 



(iV 2 > /=0 l 2 ■ P (l) + 2 2 ■ P (2) + 4 2 • P (4) = ^- + j + ^ , 
(N±)i=o = l 2 ■ -P±(l) + 2 2 -P ± (2) = ^ + t. + £!. 
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10 



(43) 
(44) 



From Eqs. ( 14111441) one finds: 
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2 
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30 



(45) 



There are no correlations between the numbers of it , tt + , tt~ in the GCE. All correlation 
coefficients defined by Eq. f l34|) are equal to zero, p + — Po- = P+~ — 0. The charge is exactly 
conserved in the Q = and 1 = statistical ensembles. This brings the strongest correlations 
between the numbers of ir + and tt~ , i.e. p + _ = 1, and this means equal numbers iV + and AL in 
each microscopic state of the system. The correlations between the numbers of 7r° and n , are 
absent in the CE with Q = 0, but exist in the statistical ensemble with 1 = 0. The correlation 
coefficient p ! 0± at I = is presented in Fig. |3l It has the maximal value Pq ± ~ 0.19 at z ~ 1 
and goes to zero at z — > oo. 
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FIG. 3: The correlation coefficient /?q ± ()34j) between the neutral and charged pions in the statistical 
ensemble with 1 = 0. 
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The role of the pion number correlations can be illustrated returning to the scaled variance 
cujv presented in Fig. HJ right. Taking into account the relation N = Nq + N + + i\L one finds: 

WK = W0 M + 2w± <M + 2/)+ _ w± i^ + 2po±v ^^. (46) 

The GCE corresponds to p + _ = p ± = in Eq. ( )46l) . Besides, the multiplicities (Nj) = (N)/3 
and scaled variances Uj = 1 are the same for j = 0, +, — . This leads to uj^ = 1 in the GCE. 
The CE with Q = corresponds to = 1, p+_ = 1 and p®± = in Eq. (j4"B"j) . This gives, 

W Q = fMM + Au Q ^)q=° s (i_ 2 ^ + 4z 2 ) + A(z - 2z 2 ) = 1 + 2z - Az 2 . (47) 

The term (N ) q =0 / (N) q =0 in Eq. (j!7|) decreases at small z. The maximum of appears due 
to pj r _ = 1, note the factor 4 in the second term in r.h.s. of Eq. (I47p . Due to the exact charge 
conservation, negative and positive pions may appear only as 7r + 7r + -pairs. For N c h = N + + A_, 
the scaled variance 

q _ ( N ch)Q=o - (Nch)Q=0 _ Q _ Q ( 

is then two times larger than that for positive (negative) pions. The same relation, u c h = 2u±, 
is also valid at 1 = 0. 

For 1 = one obtains from Eq. (146]) : 



2 z 7z 2 \ / 4 2z 2 \ /2z 11^ 2 \ 2 z 2 
3-6 + l80j + U" 45 J + U" «-J ^ 2+ 2 - 4- (49) 

where the relation (N )i =0 = (N±)i =0 = (N)i =0 /3 is used. The sum of Wq/3 and 4uj±/3 
decreases at small z. An increase of uj^ and its maximum at small z comes due to the third 



term in the r.h.s. of Eq. ( 1491) . i.e. due to the correlations between the neutral and charged 
pions, Pq ± > 0. Note that a conclusive comparison with experimental data for small pion 
systems would require an inclusion of additional conservation laws like 4-momentum, angular 



momentum, etc. 
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IV. BOSE STATISTIC 



In this section we discuss the role of Bose effects in the statistical ensemble with fixed 
isopsin. The role of Bose effects for the fluctuations in the CE was considered in Ref. 14]. A 
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generalization of Eq. fl6]) for Bose statistics gives the following expression: 



ryBose 
^1=0 z 

x exp 
where 



JO 



da dj d(3sm(3 



(50) 



o 



J2^^cos n f3 + (^ 1+ ^ OS/3 ^ (A"exp[m(a + 7 )] + A" exp[-m(a + 7)]) 



V Tm 2 



/nm\ 



(51) 



2tt 2 n " \ T / 

The Boltzmann approximation corresponds to the first term n = 1 in the sum in Eq. floPl) . The 
results for uj 1 , are presented in Fig. HI 




1 2 3 4 5 

z 

FIG. 4: The scaled variances ojq and in the Bose gas with 1 = 0. The solid lines correspond to 
m/T = 1, dashed-dotted lines to m/T -> 0, and dashed lines to m/T -> 00. 

In the case of quantum statistics, the partition function depends not only on the one particle 
partition function z 02]), but additionally on the value of m/T. In Fig. H] the results for the 
scaled variances Uj are shown for a typical value m/T = 1 and for two limiting cases: m/T — > 
when the Bose effects are the strongest ones, and m/T — )■ 00 when the Bose effects disappear 
and results coincide with the Boltzmann approximation presented in Fig. |2j right. As seen 
from Fig. HI Bose statistics makes the pion number fluctuations larger. These effects are always 
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stronger for smaller values of the m/T ratio. For m/T — > we find Uq = 2.26 at z — > 0, and 
= 1.368 at z — > oo. The corresponding results for the charged particles are the following: 
u;± = 1 at z — > 0, and bj\_ = 0.684 at z — > oo. The results for z — > oo coincide with those in the 
canonical ensemble with Q = 0. 

The Bose statistics effects become visible when more than one identical pion (i.e. pions in 
the same charge states) can appear in the same quantum state. Note a difference between the 
neutral and charged pions for very small systems with 1 = 0. At z « 1 the main configurations 
consist of 1 (negative) positive pion and 2 neutral pions. These states correspond to the total 
number of pions N = 2 and their probabilities are proportional to z 2 . The states with only 1 
neutral pion have smaller probability, z 3 , as they can only appear for N > 3. Because of this 
difference one observes no Bose effects for u± at z — > in Fig. HJ whereas these effects are seen 
for Uq. 

V. SUMMARY 

Particle number fluctuations and correlations in the statistical system of pions with zero 
total isospin 1 = have been studied in the present paper. For finite systems one observes a 
suppression of the average total pion number and an increase of the pion number fluctuations. 
The suppression effects due to isospin conservation are the same for average numbers of ir°, 7r + 
and 7r~. However, we find quite different behavior of the corresponding scaled variances. For 
neutral pions there is the enhancement of the fluctuations, whereas for charged pions the isospin 
conservation suppresses fluctuations similar to that in the canonical ensemble with Q = 0. The 
positive correlations between the numbers of neutral and (negative) positive pions are observed 
for 1 = 0. This effects is absent in the canonical ensemble with Q = 0. The correlations 
between the numbers of neutral and charged pions are responsible for a maximum of for 
small systems with 1 = 0. In the thermodynamic limit the correlation coefficient p T Q± between 
neutral and charged pions goes to zero. This leads to Wq — > 1, u± — > 1/2 for Boltzmann 
statistics and to Uq —> 1.368, <jj± — » 0Jq/2 = 0.684 for m = and Bose statistics. 

Finally, we speculate on possible consequences of our results for the quark-gluon gas with 
SU(3)-color group. A requirement of colorless of the system of quarks and gluons is similar 
to that of the isospin singlet 1 = 0. This requirement cause the suppression effects for the 
average yields of quarks and gluons, the same for different colors. These suppression effects 
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for the yields disappear in the thermodynamic limit. However, a behavior of the fluctuations 
are different. To make situation fully transparent let us consider a toy model of the gluon gas 
with SU(2)-color. The colorless system of gluons will be then in one-to-one correspondence 
with the iso-singlet 1 = pion gas. We thus immediately conclude that fluctuations of the 
number of gluons with different colors are different. Moreover, this difference survives in the 
thermodynamic limit too. 
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